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Abstract. We study Maxwell's equations on a 4-inanifold where the elec- 
tromagnetic medium is modelled by an aiitisyniinctric (^j-tensor with real 
coefficients. In this setting the Presnel surface is a fourth order polynomial 
surface in each cotangent space that acts as a generalisation of the light cone 
determined by a Lorentz metric; the Presnel surface parameterises electromag- 
netic wave-speeds as a function of direction. The contribution of this paper is 
a pointwise description of all electromagnetic medium tensors that satisfy the 
following conditions: 
(i) K is invcrtiblc, 
(ii) K is skcwon-free, 

(Hi) K is birefringent, that is, the Presnel surfa<;e of K is the union of two 
distinct light cones. 

Wc show that there arc only three classes of mediums with these properties 
and give explicit expressions in local coordinates for each class. 



We will study the pre-metric Maxwell's equations. In this setting Maxwell's equa- 
tions are written on a 4-manifold and the electromagnetic medium is described 
by an antisymmetric (2) -tensor k on N. Then the electromagnetic medium k 
determines a fourth order polynomial surface in each cotangent space called the 
Fresnel surface ^ and it acts as a generalisation of the light cone determined by 
a Lorentz metric; the Fresnel surface parameterises wave-speeds as a function of 
direction [Rub02, HO03, PSW09]. At each point in spacetime -/V, the electromag- 
netic medium depends, in general, on 36 free components, In this work we assume 
that the medium is skewon-free. Then there are only 21 free components and such 
medium describe non-dissipative medium. For example, in skewon-free medium 
Poynting's theorem holds under suitable assumptions. 

The above means that in the pre-metric setting we have two descriptions of elec- 
tromagnetic medium: First, we have the (2) -tensor k that contains the coefficients 
in Maxwell's equations. On the other hand, we also have the Frcisnel surface 
which describes the behaviour of a wavespeed for a propagating electromagnetic 
wave. If K is known we can always compute ^ by an explicit equation (see equa- 
tion (10)). A more challenging question is to understand the converse dependence, 
or inverse problem: If the Fresnel surface is known for some p & N, what can 
we say about k|p? Essentially this asks that if the behaviour of wave speed for 
an electromagnetic medium is known, what can we say about the medium? These 
questions are of theoretical interest, but also of practical interest as they relate to 
understanding measured data in engineering applications like traveltime tomogra- 
phy in anisotropic medium. We will here only study the problem at a point p & N 
since the dependence will never be unique. For example, the Fresnel surface ^ is 
always invariant under scalings and inversions of k [HO03, Dahlia]. In general, 
these are not the only invariances, and for a general k the relation between k and 
does not seem to be very well understood. 
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A natural first task is to characterise those mediums k for which the Fresnel surface 
^ is the light cone of a Lorentz metric g. This question was raised in [OH99, 
OFROO]. A partial solution was given in [OFROO], and (in skewon-free mediums 
with real coefficients) the complete solution was given in [FBll]. The result is that 
if the Fresnel surface is a light cone, then k is necessarily proportional to a Hodge 
star operator (plus, possibly, an axion component proportional to the identity). For 
an alternative proof, sec [Dahlia], and for related results, sec [OR02, LH04, Iti05]. 

The main contribution of this paper is Theorem 2.1. It gives a pointwise charac- 
terisation of all electromagnetic medium tensors k with real coefficients such that 

(i) K is invertible, 

(a) kSs skewon-frcc, 

(in) K is birefringent, that is, the Fresnel surface of k is the union of two distinct 
light cones for Lorentz metrics. 

The first two assumptions imply that k is essentially in one-to-one correspondence 
with an area metric. Area metrics also appear when studying the propagation of 
a photon in a vacuum with a first order correction from quantum electrodynam- 
ics [DH80, SWWIO]. The Einstein field equations have also been generalised into 
equations where the unknown field is an area metric [PSW07]. For further exam- 
ples, see [PSW09, SWWIO] and for the differential geometry of area metrics, see 
[SW06, PSW07]. For Maxwell's equations, the interpretation of condition (Hi) is 
that diflFerently polarised waves can propagate with different wave speeds. In such 
medium one should expect that propagation of electromagnetic waves is determined 
by null-geodesics of two Lorentz metrics. A typical example of such medium is a 
uniaxial crystal For partial results describing when the Fresnel surface factorises, 
see [RSll] and in 3 dimensions, see [Kac04, Dah06]. 

In Theorem 2.1 wo show that there arc only three medium classes with the above 
properties and we give explicit expressions in local coordinates for each class. Of 
these classes, the first is a generalisation of uniaxial medium and the last seems 
to be unphysical; heuristic arguments suggest that Maxwell's equations are not 

hyperbolic in the last class. 

The main idea of the proof is as follows. We will use the normal form theorem 
for area metrics derived in [SWWIO], which pointwise divides area metrics into 23 
metaclasscs and gives explicit expressions in local coordinates for each metaclass. 
This result was also used in [FBll], and by [SWWIO] we only need to consider the 
first 7 metaclasses. For each of these metaclasses, the Fresnel surface can be written 
as = {(, EM.^ : f{£,) = 0} for a homogeneous 4th order polynomial / : — >■ K 
with coefficients determined by K\p. Since k\p is birefringent, / factorises as 

/(O = ^eM^ 

into homogeneous 2nd order polynomials /± : — >■ M. By identifying coefficients 
we obtain a system of polynomial equations in coefficients of / and f±. In the 
last step we eliminate the coefficients in f± from these equations whence we obtain 
constraints on / (and hence on k) that much be satisfied when k is birefringent. 
To eliminate variables we use the technique of Grobner bases, which was also used 
in [Dahlia]. 

A limitation of Theorem 2.1 is that the explicit expression is only valid at a point. 
The reason for this is that the decomposition in [SWWIO] essentially relies on the 
Jordan normal form theorem for matrices, which is unstable under perturbations. 
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Another limitation is that we do not aUow for complex coefficients in n. Therefore 
mediums like chiral medium are not included in the mediums in Theorem 2.1. 

This paper relies on computations by computer algebra. For information about the 
Mathematica notebooks for these computations, please see the author's homepage. 

1. Maxwell's equations 

By a manifold M we mean a second countable topological Hausdorff space that 
is locally homeomorphic to M" with C°°-smooth transition maps. All objects are 
assumed to be smooth and real where defined. Let TM and T* M be the tangent 
and cotangent bundles, respectively. For /c > 1, let A'^(M) be the set of antisym- 
metric fc-covectors, so that k^{N) = T*N. Also, let nf(M) be (;)-tensors that 
are antisymmetric in their k upper indices and I lower indices. In particular, let 
n''{M) be the set of /c-forms. Let C°°{M) be the set of functions. The Einstein 
summing convention is used throughout. When writing tensors in local coordinates 
we assume that the components satisfy the same symmetries as the tensor. 

1.1. Maxwell's equations on a 4- manifold. On a 4-manifold N, Maxwell's 
equations read 

(1) dF = 0, 

(2) dG = j, 

where d is the exterior derivative on N, F,G G ^"^{N) are the electromagnetic field 
variables, and j G fl^{N) is the source term. By an electromagnetic medium on N 
we mean a map 

We then say that 2-forms F,G & Q'^{N) solve Maxwell's equations in medium k if 
F and G satisfy equations (l)-(2) and 

(3) G = >.{F). 

Equation (3) is known as the constitutive equation. If k is invertible, it follows 

that one can eliminate half of the free variables in Maxwell's equations (l)-(2). We 
assume that k is linear and determined pointwise so that we can represent k by an 
antisymmetric (2)-tensor k G Ct2i^)- '^^ coordinates {a;'}?^o for N we have 

(4) K = ^'^^i-^'^'®^^"® ^® 

and F = Fijdx^ (S) dx^ and G = Gijdx^ ® dx^ , then constitutive equation (3) reads 

(5) = 

Then at each point on N, a general antisymmetric (2) "tensor k depends on 36 free 
real components. In the main result of this paper (Theorem 2.1) we will assume 
that K is skewon-free, that is, k{u) Av = uA k{v) for all u,v G 0,2{N) whence 
K has only 21 free components. Physically, such medium describe non-dissipative 
medium; if k is time-independent, then Poynting's theorem holds under suitable 
assumptions [DahlO, Proposition 3.3]. Let us also note that if A'' is orientable, then 
invertible skewon-free mediums are essentially in a one-to-one correspondence with 
area metrics. See [SWWIO] and [Dahllb, Proposition 2.4]. The medium is called 
axion-free if trace k = [HO03] . 

By a pseudo-Riemann metric on a manifold M we mean a symmetric (2) -tensor g 
that is non-degenerate. If M is not connected we also assume that g has constant 
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signature. By a Lorentz metric we mean a pseudo-Riemann metric on a 4-manifold 

with signature (— + ++) or (H ). The light cone of a Lorentz metric is defined 

as 

N{g) = {CgT;(7V): 5(^,0=0}. 
For pe N wc define Np{g) = N{g) T*{N). 

ff 5 is a pseudo-Riemann metric on a orientable 4-manifold N, then the Hodge 
operator of g induces a skewon-free (j) -tensor that we denote by k = *g. Moreover, 
if locally g = gijdx^ dx\ and k is written as in equation (4), then 

(6) <i = V\9\9"'9'''eabrs, 

where detg — dot g-y, g*^ is the ijth entry of {gij)~^, and £(i...;„ is the Levi-Civita 
permutation symbol. We treat Si^-.-i^ as a purely combinatorial object (and not as 
a tensor density). Let also g'l'"'* = ei^...i^. 

1.2. Representing k as a 6 x 6 matrix. Let O be the ordered set of index pairs 
{01, 02, 03, 23, 31, 12}. If 7 e O, we denote the individual indices by h and h- Say, 
if 7 = 31 then h = 1. 

If {x'}^^ Q are local coordinates for a 4-manifold TV, and J G O we define dx = 
dx-^^ A dx-^^. A basis for n'^{N) is given by {dx-^ : J G O}, that is, 

(7) {dx° A dx^,dx° A dx"^, dx° A dx^, dx^ A dx^, dx^ A dx^.dx^ A dx^}. 

This choice of basis follows [HO03, Section A.1.10] and [FBll]. If k e ^1{N) is 
written as in equation (4) and J G O, then 

(8) K{dx-^) = ^ Kjdx^, J eO, 

where jj- We will always use capital letters I, J,K,. . . to denote elements 

in O. Let b be the natural bijection 6: O — >■ {1, . . . , 6}. Then we identify coefficients 
{k/ : 7, J e O} for k with the 6x6 matrix A = {Kj)ij defined as nj = 
for 7, Je O [Dahllb]. 

1.3. Presnel surface. Let k <E f^ll^) ^ 4-manifold N. If k is locally given by 
equation (4) in coordinates {a;*}|^Q, let 

C^ijkl _}^^aia2^a3i atj bib2b5k b3b4b6l 

~ 48 i'il>2 '^babi'^bsbe^ ^ £01020304- 

In overlapping coordinates {x^}^^q, these coeflacients transform as 

~ dx- dx^ dx- dx'^' 

and components {^J-'*^' define a tensor density % on A'' of weight 1. The Tamm- 
Rubilar tensor density is the symmetric part of 5fo and we denote this tensor density 
by ^ [Rub02, HO03, Dahlia]. In coordinates, = where parenthesis 

indicate that indices ijkl are symmetrised with scaling 1/4!. Using tensor density 
^, the Fresnel surface at a point p G N is defined as 

(10) = U G Ai(A^) : = 0}. 

The Presnel surface is a fundamental object when studying wave propagation in 
Maxwell's equations. It is clear that in each cotangent space, the Fresnel surface 
is a fourth order polynomial surface, so it can have multiple sheets and sin- 
gular points. There are various ways to derive the Presnel surface; by studying 
a propagating weak singularity [OPROO, Rub02, IIO03], using a geometric optics 
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[Iti09, Dahlia], or as the characteristic polynomial of the full Maxwell's equations 
[SWWIO]. Classically, the Presnel surface can be seen as the dispersion equation 
for a medium, so that it constrains possible wave speed(s) as a function of direction. 

If K = f*g for a Lorentz metric g and a non-zero function / G C°°(N), then the 
Fresnel surface is the light cone of g. The converse is also true (assuming that k is 
skewon-free and axion-free) [HO03, FBll, Dahlia]. The medium given by k = f*g 
is known as non-birefringent medium. For such medium propagation speed does 
not depend on polarisation. 

Definition 1.1. Suppose is a 4-manifold, k S n2{N) and ^ is the Fresnel 
surface for k. Up G N we say that decomposes into a double light cone if there 
exists Lorentz metrics g+ and g- defined in a neighbourhood of p such that 



and Np{g+) ^ Np{g-). We then also say that k\p is birefringent. 

We know that two Lorentz metrics are conformally related if and only if they have 
the same light cones [Ehr91]. The condition Np{g+) ^ Np{g-) thus only exclude 
non-birefringent mediums, which for skewon-free mediums are well understood (see 
above). When ^\p decompose into a double light cone as in Definition 1.1 a physical 
interpretation is that the medium is birefringent. That is, differently polarised elec- 
tromagnetic waves can propagate with different wave speeds. Common examples 
of such mediums are uniaxial crystals like calcite [BW80, Section 15.3]. 

To prove of the next four propositions we will need some terminology from al- 
gebraic geometry. If A; = K or A; = C, we denote by k[xi,. . . ,Xn] the ring of 
polynomials fc" — >■ /c in variables Xi, . . . , a;„. Moreover, a non-constant polynomial 
/ e fc[xi, . . . , Xn] is irreducible \i f = uv for u,v & k[xi, . . . , Xn] implies that w or t; 
is a constant. For a polynomial r G k[xi, . . . , a;„], let V{r) = {a; S fc" : r{x) = 0} 
be the variety induced by r, and let (r) = {/r : / e k[xi, . . . ,Xn]} be the the 
ideal generated by r. For what follows the necessary theory for manipulating these 
objects can, for example, be found in [CL092]. 

The next proposition can bo seen as a reformulation of the Brill equations that 
characterise when a homogeneous polynomial factorises into linear forms [BrilOj. 



Proposition 1.2. Suppose Q G c^^** is a symmetric non-zero matrix and f G 
C[Co, • • • , 6] is the polynomial /(^ = C ■ Q ■ Uor ^ = (Co, . • . , 6) € Then f is 
irreducible in C[^o, ... , ^3] if and only if adj Q ^ 0. 

In Proposition 1.2, adj Q is the adjugate matrix of all cofactor expansions of Q, and 
^* is the matrix transpose. 

Proof. The result follows from the following three facts: First, ii f = uv for poly- 
nomials u,v G C[^o, • • • jCs] then u and v are linear. (To see this, we know that 
(5(0) = 0, so we may assume that m(0) = 0. For a contradiction, suppose that 
v{0) ^ 0. Then (i/lo = implies that du\o = 0, but then u = 0, so Q = 0.) Second, 
by [BrilO, Example 2] polynomial / is a product of linear forms or / has a multi- 
ple factor if and only if G/(^, 77, C) = for all ^,r],( G C^. Here G/ is the Gaeta 
covariant defined as 



(11) 



^\p = Np{g+) U Np{g_) 




2/(?) {Df)dv) {Df)dC) 
iDf)^{n) 2f{rj) {DfUO 

(D/)j(C) {DfUO 2/(0 
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and {Df)a{b) = ^(/(a + <6))|t=o is the directional derivative. Third, by computer 
algebra we have adj Q = if and only if G/ = 0. □ 

In [Mon07] it is proven that the light cone of a Lorentz metric can not contain a 
vector subspace of dimension > 2. The next proposition generalise this result to 
double light cones. In the proof of Theorem 2.1 we will use this result to show 
that medium tensors in the last 16 metaclasscs in the classification of [SWWIO] 
can not decompose into a double light cone. In [SWWIO] this property was used 
to show that these last mctaclasses are neither hyperbolic, so in these metaclasses, 
Maxwell's equations are not well-poscd. 

Proposition 1.3. Suppose g± are Lorentz metrics on a A-manifold N. //F C TpN 
is a non-empty vector subspace such that T c N{g^) U N{g_), then dimF < 1. 

Proof. Wc may assume that dimF > 1. Let us first prove the result in the spe- 
cial case that g+ and (/_ are conformally related (after [Mon07, Proposition 2]). 
Let {a;*}f^Q be coordinates around p such that g^\p = ±diag(— 1, 1, 1, 1), whence 
we can identify TpM and M © M"^. Let tt be the Cartesian projection onto the 
first component. Then the restriction 7r|r: F — )• R satisfies kerTrjr = {0} and 
dimrange7r|r = 1, and the claim follows. 

For general g+ and g- let us show that dimF > 2 leads to a contradiction. If 
dimF > 2, we can find linearly independent u,v £ T such that spa.n{u,v} C 
N{g^) U N{g-). We may further assume that u G N{g+). Let 

U = {6 €R: coseu + sinev ^ N{g_)}. 

For w G T* N let us write = g+iw, w). If U is empty, then spanjii, v} C N{g^) 

and the result follows from the special case. Otherwise there exists a G so 
that II cosOu + sinOvW^ = for all 9 in some neighbourhood Iq 3 Oq. Differentiating 
gives 

^ {\\v\\^ -\\u\\^) ■ sm2e + g+{u,v) ■ cos2e = 0, 6 & Iq. 

By computing the Wronskian of sin 20, cos 2^, it follows that = ||u||^ = ||u||^ and 
g-^.{u,v) = 0. Thus span{u, c N{g+), but this contradicts the special case. □ 

The next proposition gives the pointwise description of the Tamm-Rubilar tensor 

density at points p E N where the the Fresnel surface decomposes into a double 
light cone. Let us emphasize that the result is pointwise. For example, in equation 
(12) the two sides have different transformation rules. 

Proposition 1.4. Suppose N is a A-manifold, k G i^i^); '^^'^ ^he Fresnel surface 

of K decomposes into a double light cone at p G N . If {s'jf^g are coordinates 
around p and '^'^^^^ and g± = g±ijdx^ (E) dx^ are as in Definition 1.1, then 

(12) '^'^'''UMi = C{g%i,){g''J^k^i) atp, {^j^o e M^ 

for some C € M\{0}. 

Proof Let /±,7: ^ C be polynomials /±(0 = gtUi and 7(0 = ^^^'^'6066 
for ^ = (^i)f^Q e C^. For ideals I± = (/±) we then have {f+f-) = I+H I- whence 
equation (11) implies that V{{'y)) = V{I+ fl /_) and passing to ideals gives 

mm = i{v{uni.)). 

The Strong Nullstellensatz implies that [CL092, p. 175] 

(7) C ^/^n^/^, 
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where \fl be the radical of an ideal / and we used identities a/iTTJ = Vl n VJ 
and I C Vl valid for any ideals I and J [CL092, Proposition 16 in Section 4.3]. 
By Proposition 1.2, f± are irreducible polynomials, so I± are prime ideals and 
I± = \fl±- Thus 

(7> C (/+/_) 

and 7 = pf+f- for some polynomial p: — > C. Computing the polynomial degree 
of both sides implies that p is a non-zero constant p = C € C. By Proposition 1.3, 
we can find a ^ e R^* so that ^ ^ V{f+) U V{f_) = V{'y) whence C € K\{0}. □ 

For two Lorentz metrics g and h we know that their light cones N{g) and N{h) 
coincide if and only if g and h are conformally related [Ehr91]. The next proposition 
gives an analogous uniqueness result for Fresnel surfaces that decompose into a 
double light cone. 

Proposition 1.5. Suppose N is a 4-manifold, k G 0,2{N) and the Fresnel surface 

of K decomposes into a double light cone at p € N . Suppose furthermore that 
equation (11) holds for Lorentz metrics g± and for Lorentz metrics h±. Then there 
exists constants C± G 1R\{0} such that exactly one of the following conditions hold: 
g± = C± h± at p or gzf = C± h± at p. 

Proof. The result follows by Propositions 1.2 and 1.4 and since any polynomial 
has a unique decomposition into irreducible factors [CL092, Theorem 5 in Section 
3.5]. □ 

The next example shows that unique decomposition is not true when ^ only de- 
compose into second order surfaces. 

Example 1.6. In coordinates {x^}f^Q for R^, let k be the skcwon-free (jj'tensor 
determined by the 6x6 matrix = diag(— 1, 1, 0, —1, 1, 0). Then k has Fresnel 

surface 

^ = e t*k4 : ^oCi66 = 0}. 

It is clear that ,^ has multiple factorisations into quadratic forms, and by Propo- 
sition 1.3, ^ does not factorise into a double light cone. □ 

2. NON-DISSIPATIVE MEDIA WITH A DOUBLE LIGHT CONE 

Theorem 2.1 below is the main result of this paper. To formulate the theorem we 
first need some terminology. Suppose L: V ^ V ]s & linear map where V is & n- 
dimcnsional real vector space. If the matrix representation of L in some basis is A G 
j^nxra g^j^j -^yj-j^i^eji using the Jordan normal form we say that L has Segre type 
[toi • • • TOr k\k\ - ■■ ksk^ when the blocks corresponding to real eigenvalues have 
dimensions mi < • • • < and the blocks corresponding to complex eigenvalues 
have dimensions 2ki < ■ ■ ■ < 2kg. Moreover, by uniqueness of the Jordan normal 
form, the Segre type depends only on L and not on the basis. For a (2) -tensor 
K on a 4-manifold, we define the Segre type of k\p as the Segre type of the linear 
map n'^{N)\p — n'^{N)\p in basis (7). By counting how many ways a 6 x 6 matrix 
can be decomposed into Jordan normal forms, it follows that there are only 23 
Segre types for a (2) -tensor. A main result of [SWWIO] arc simple normal forms 
in local coordinates for each of these Segre types under the assumption that k\p is 
skewon-free and invertible. In the below proof we will use the restatement of this 
result in [Dahllb]. 
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In the proof of Theorem 2.1 we will eliminate variables in systems of polynomial 
equations. Suppose F C C" is the solution set to polynomial equations /i = 
0, . . . , /jv = where fi G C[xi, . . . , a;„]. If I is the ideal generated by /i, . . . , /at, 
the elimination ideals are the polynomial ideals defined as 

Ik = /nC[a;fe+i,...,a;„], fc G {0, . . . , n - 1}. 

Thus, if (xi, . . . , Xn) € V then p{xk+i, . . . , Xn) = for any p € Ik, and Ik contain 
polynomial consequences of the original equations that only depend on variables 
Xh+i, . . . ,Xn- Using Grobner basis, one can explicitly compute Ik [CL092, Theorem 
2 in Section 3.1]. In the proof of Theorem 2.1, we will use the built-in Mathematica 
routine GroebnerBasis for such computations. 

Theorem 2.1. Suppose N is a 4-manifold and k S f2|(A^). Furthermore, suppose 
that at some p G N 

(i) K,\p has no skewon component, 
(a) k\p is invertible as a linear map A.p{N) Ap{N), 
(Hi) the Fresnel surface factorises into a double light cone. 

Then k\p must have Segre type [iT iTlI], [22 iT] or [11 iTlT]. 

(i) Metaclass I: If k\p has Segre type [11 11 iT], there are coordinates {a;*}?^o 
around p such that 
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(13) A = 

where i' and i" are defined such that {i,i',i"} = {1,2,3} and i' < i" . 
Then, for exactly one i G {1, 2, 3} we have 

(14) = 2, A' = A" 
and equation (11) holds for Lorentz metrics 



9± 



diag I 1 
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-D,±JDi-4 , 
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-A ± a/A -4 



-A3 ± ^ Di - 



(a) Metaclass II: If k\p has Segre type [22 11], If k\p is in Metaclass II, there 
are coordinates {a;*}?^g around p such that 
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where ai,a2 S K and Pi > 0. Then o.\ = 0.1 and /3i = l32, and equation 
(11) holds for Lorentz metrics 



9± = 



( ±1 








/3i 







-Pi 
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(Hi) Metaclass IV: If k\p has Segre type [111111], there are coordinates 
{a;'}f^o around p such that 



( Oil 








-/3i 










0:2 








-P2 











013 








ai 


Pi 








ai 











P2 








0:2 





^0 





ai 








013) 



for some ai, Q!2, as, a4 € IR anrf /3i,,52 > 0. Then ai = a^, Pi = P2, 
ai ^ 0, 03 ^ a\ and equation (11) holds for Lorentz metrics 

9± = diag (^1, ^ ±^Dl+4j,^ (^^Di ± ^Dl + 4^ , -1^ , 

w/iere 

_ (aa - a3)2 + ^2 _ ^2 



Proof. Metaclass I. The local expression for k\p is given by [Dahlia, Theorem 
3.2]. Then the Tamm-Rubilar tensor density for k\p satisfies 

(16) C-^'^'i'^'^i^i^k^i = ^o^ + e^ + d+Cs*- ^0^0^166 

3 

i=l 

where C = P1P2P3 and Dq is given explicitly in terms of ai, . . . , ^3, and implicitly 

Df) satisfies 

(17) Df,=A{4 + DiD2D:,-Dl-Dl-Dl). 

By Proposition 1.4, there are real symmetric matrices A = {A^^)^j_q and B = 
{B'^)lj^o such that 

(18) C-'^r^'^'^jCk^i = (e-A-^) ie-B-^), ^eR\ 

Writing out these equations shows that A^^S*"^ = 1. Hence is non-zero, and 
by rescaling A and B, we may assume that A^^ = 1. This substitution simplifies 
the equations so that by polynomial substitutions wc may eliminate all variables 
in B and variable Dq. This results in a system of polynomial equations that only 
involve Di,D2,D^ and the variables in A. Further eliminating the variables in A 
using a Grobncr basis, gives constraints on Di. D2, D:^. By equation (13) wc know 
that Di, D2, D3 > 2, whence these constraints imply that there exists a unique 
i e {1, 2, 3} such that conditions (14) hold. (To see that i is unique it suffices to note 
that if Di = Dj = 2 for i 7^ j then Di = D2 = D3 = 2 whence Dq = and equation 
(18) holds for the single light cone A = B = g^^ with go = diag{— 1, 1, 1, 1}. 
By Proposition 1.4 and unique decomposition into irreducible factors this gives a 
contradiction.) Equations (14) and (17) imply that Dq = 0. The result follows 
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since equation (18) holds with A = and B = g_ , where g± are the matrices 
in the theorem formulation. 

Metaclass II. As in Metaclass I, there are matrices A and B such that the Presnel 
polynomial satisfies equation (18) (with C = 1). As in Metaclass I we can eliminate 
variables in B. Further eliminating all variables in ^ by a Grobner basis implies 
that ai = a2 and /?i = (32- Then a a direct computation shows that equation (18) 
holds with A = g^^, B = gZ^ and C = j3\. Computer algebra shows that g± have 
Lorentz signatures. 

Metaclass III. If k\p is in Metaclass III, there are coordinates {a;*}?^g around p 
such that 



(oil 


-A 











\ 


Pi 
















1 





ai 








-Pi 











ai 


Pi 


1 








1 


-Pi 


ai 





^0 


1 


Pi 








ai ) 



where ai S M and j3i > 0. Decomposing the Frcsncl polynomial as in equation (18) 
(with C = 1) gives a system of polynomial equations for the variables in A, B and 
K\p. Computing the Grobner basis for these equations implies that P\ = 0. Thus 
k\p can not be in Metaclass III. 

Metaclass IV. Let us first note that 0:4 7^ since otherwise span{dx^\p,dx'^\p} C 
^\p which is not possible by Proposition 1.3. Then the Tamm-Rubilar tensor 
density satisfies 

+^1(^1^3' - eofi) + D^ieies - eoe2) + i?3(-e??i - ^o^^^), 

where C = P1P2O.4, Dq is determined explicitly in terms of ai, . . . , /32, Di is defined 
in equation (15), IPs > 2 is defined in equation (13) and 

(cti - as)' + Pf - aj 

L>2 = a • 

PlQ!4 

By decomposing and eliminating variables as in Metaclass I, it follows that that 

Dq = and D-^ = 2. Thus wc have proven that ai = 0.2 a-nd /3i = /32 whence D\ = 
D2 and equation (18) holds with A = gZ^ , B = gZ} and C as above. Moreover, g± 
both have Lorentz signatures. Condition ^ a\ follows since det k\p ^ 0. 

Metaclass V. If k\p is in Metaclass V, there are coordinates {a;'}f^o around p such 
that 



(oii 


-Pi 











o\ 


Pi 






















a2 








as 





1 





ai 


Pi 





1 








-Pi 


Oil 





^0 





az 








Q!2/ 



where ai,a2,as G M and /3i > 0. We may assume that as ^ 0, since otherwise 
span{da;*|p}?^i C ^\p- Decomposing and eliminating variables as in Metaclass I 
gives that Pi is purely complex. Thus /t|j, can not be in Metaclass V. 
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Metaclass VI. If k\p is in Metaclass VI, there are coordinates {a;*}|^o around p 
such that 











-01 










a2 








a4 











as 










/3i 








a I 











ai 








a2 





\0 





Q!5 








"3/ 



/J 



for some ai , . . . , as e M and > 0. We may assume that a4 and are non-zero 

since otherwise spanjda;*, dx'^} C J^\p for some i € {0, 1} as in Metaclass IV. Then 
the Tamm-Rubilar tensor density satisfies 

where C = Pia^a^ and Dq, Di, I?2i G M are defined in terms of a^ and By 
decomposing the Fresncl tensor as in equation (18) and eliminating variables using 
a Grobner basis, it follows that there exists a a € {±1} such that 

Do = 0, Di= cr2, Da = -(7D3, 

and moreover, equation (18) holds ioi A = g^^ , B = gZ^ and C as above, where 

1 



g± = diag l,-f7, - <t£>3± 



Since g+ does not have a Lorentz signature for any a e {±1} and Da G M, Propo- 
sition 1.4 and unique factorisation imply that K\p can not be in Metaclass VI. 

Metaclass VII. If k\p is in Metaclass VII, there are coordinates {a;'}f^o around p 
such that 



/ ai 








a4 





o\ 





a2 








as 











as 








ttfj 


a4 








ai 











0:5 








a2 





^0 





ae 








Q!3/ 



/J 



for some ai,...,a6 S ffi. We may assume that a4,a5,a6 7^ since otherwise 
spa,n{dx^\p,dx^\p} C for some i,j e {0,1,2} as in Metaclass IV. Then the 
Tamm-Rubilar tensor density satisfies 

3 

C-l = ^0' + ^1 + el + ^3 + ^0C0666 - E ^^(e''e'" + ^0&), 



i=l 



where C = a^^a^aQ, constants Di, D2, D3 e M arc given by 

(a2 - as)^ - Oil- al 



(19) 
(20) 
(21) 



(ai 



as 



asae 

2 



al 



afi 



a4af; 
(ai - aa)^ - 



9 

a| 



a4a5 



and Do e M is given explicitly in terms of ai, . . . , /^s, and implicitly Dq satisfies 
(22) Dl=4{-4 + DiD2D3 + Dl + D^ + Dl) . 
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Decomposing the Tamm-Rubilar tensor density as in equation (18) and eliminating 
variables using a Gr5bner basis, gives polynomial equations for Do,Di,D2,Ds. Let 
us consider the cases Dq = and Dq separately. If Dq = 0, there exists an 
i e {1, 2, 3} and a a G {±1} such that 



(23) 



£>o = 0, Di = -a2, Di,=aDi„, 



where the last condition is a consequence of equation (22). Suppose i = 1. Then 
Proposition 1.4 implies that for some invertible symmetric matrices A,B€ R*^^ 
with Lorentz signatures we have 

(24) {e-A-^){e-B-c) = {e-L+-^){e-L_-^), ^ec^ 

where matrices L± g C^^'* are defined as 

(25) L± = diag(^l,a,i(^-A±/Di^),|(^-i}2±/Dr^)). 

Since L± are invertible, equation (24), Proposition 1.2 and unique factorisation 
imply that L± are real and have Lorentz signatures. Thus \D2\ > 2 and det L± < 0, 
but this contradicts equation (25), which implies that 



det 



L± = ^(-D2±^Dl-4^ >0. 



A similar analysis for i = 2, 3 shows that the case JDq = is not possible. If Dq 7^ 
it follows that there exists ai, 0-2, G {±1} and distinct k G {1, 2, 3} such that 

(26) Do ^0, A = CTi2, Dj = £722, -Dfe = ^ (-4£7ia2 + asA) , 

where the last equation follows from equation (22). If = (1, 2) then fc = 3 and 
Proposition 1.4 implies that for some invertible symmetric matrices A,Bg R*^^ 
with Lorentz signatures, equation (24) holds for matrices L± G C*^^ defined as 



(27) 



L+ = 



-£7l 







Vs 

-(72 





± ^ \ 







(71(72 



Since both sides in equation (24) should decompose into the same number of ir- 
reducible factors, it follows that ^* • L± ■ ^ are irreducible in C[^07 • • • 7 Cs]- Thus 
equation (24) and unique factorisation imply that L± are real and have Lorentz 
signatures, so det L± < 0. However, this contradicts equation (27) which implies 
that 



detL± = Qa - cria2cr3^ > 0. 
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The cases = (1, 3), (2, 3) are excluded by the same argument by using metrics 



/ 



L± 



1 






/ 1 



L± 



± 





-CTl 



/D7, 



±- 



± 



V 








CTiCr2 






V8 



C7lCr2 





-at 

V8 









-0-2 




-0-2 J 



respectively. Thus k\p can not be in metaclasses VII. 

Metaclasses VIII— XXIII. (Following [SWWIO, Lemma 5.1].) Let A = 
be the 6x6 matrix that represents k\p in some coordinates {a;*}f^g around p. 
Then the Jordan normal form of A has a block of dimension d € {2, . . . , 6} that 
corresponds to a real eigenvalue A G M\{0}. By considering unit vectors in the 
normal basis, we can find non-zero 61,62 S A^(A'')|p so that K{e\) = Aei and 
k(62) = Ae2 + ei. Writing out K(ei) A 62 = ei A «;(e2) implies that ei A ei = 0, so 
61 = ryi A 772 for some linearly independent r?i,7?2 € A^(iV)|p [Coh05, p. 184]. Let 
W = span{?7i, 772}- For all ^ e W we then have 

whence Theorem 3.3 in [Dahlia] implies that W C ^\p and Proposition 1.3 implies 
that k\p can not be in metaclasses VIII- XXIII. □ 



In the proof of Theorem 2.1 the assumption that k is invertible is only used to 
show that ^ a1 in Metaclass IV and to exclude Metaclasses VIII — XXIII. It 
would be interesting to sec if these last metaclasses can be excluded also for non- 
invertible k by other arguments. Regarding this question it should be emphasized 
that here k is real. For complex coefficients k, the setting becomes more involved. 
For example, in [Dahlia, Example 5.3] it is shows that for complex k the Fresnel 
surface can be a single light cone even if k is not invertible. Also, chiral medium 
would be an example of a medium with complex coefficients in k and with a double 
light cone. In chiral medium right and left hand circularly polarised plane waves 
propagate with different wavespeeds. Let us note that if we set ^0 = in equation 
(16) we obtain the ternary quartic studied in [Thol6] and for this polynomial, Dq 
in equation (17) is one of the factors in the discriminant. 

Let us make some comments regarding the three mediums derived in Theorem 2.1. 
A first observation is that for each Metaclass in Theorem 2.1, the light cones are 
parameterised by only one parameter: Dii = Di" > 2 in Metaclass I, /3i > in 
Metaclass II, and Di € K in Metaclass IV. Let consider each metaclass under the 
assumption that Theorem 2.1 holds. 

Metaclass I. In Metaclass I we assume that conditions (14) holds for only one 
i e {1, 2, 3}. In terms of ai, . . . , /^s conditions (14) are equivalent to 

When ai = a2 = cts = this medium reduces to a uniaxial medium, where wave 
propagation is well understood. 
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Let US also note that if Di = Dj = 2 for two distinct i,j E {1,2,3}, then ai ~ a2 ^ 
Q3 and Pi ~ (^2 — h- Then k|p = — /3i *g +ai Id for the locally defined Lorentz 
metric g = diag(— 1, 1, 1, 1) and the Fresnel surface is the single light cone 

Metaclass II. For Metaclass II, the Fresnel polynomial ^^^-''^'CiCjCfcCi is a function 
of ?OiCi +^2t^3- It is therefore motivated to project onto = 0, and we can 
plot as a surface in M.^ . Figure 1 show this projection for three different values 
of /3i. From the figures (or from metrics g±) we see that the light cones N{g±) 
coincide in the limit /?i — > cxd. 



Figure 1 . Projection into of Fresnel surfaces in Metaclass II 
for /3i 0.2 (left), /3i = 0.8 and /3i = 5 (right). 



By a coordinate transformation we can put the local representation of k\p into a 
more symmetric form. Let {P}f^o coordinates defined as P = J2^=o ^ij^^ 
where L — (Lij) is the Jacobian matrix 

/O ^^{l-w) 4(l + ^)\-^ 
.^01 
10 
\0 1 1 / 

where w — \/l + 4/3j. The motivation for these coordinates is that they diagonalize 
Then the 6x6 matrix (iij)ij that represents n\p in {Pjf^o coordinates is 



Qi Id H — 



/ 











PI 








\ 







Pi 


-Pi 





Pi{-l + w) 


-Pi 









Pi 


-Pi 





-Pi 


-Pl{l + 


w) 




-w 
























~Pi{l + w) 


Pi 





Pi 


Pi 




I 





Pi 


Pi{-l + w) 





-Pi 


-Pi 


) 



2.1. Metaclass IV. For Metaclass IV, the Fresnel polynomial is also a function of 
^o,£,i + £.2 1^3- We may therefore visualize the Fresnel surface in the same way as 
in Metaclass II. See Figure 2. 
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Suppose tti = a2 = ^3 = 0. If we treat x° as time and {x^}i^i as space coordinates, 
and write Maxwell's equations using vector fields E, D,B,H, then k\p represents 
medium 

D -diag(;3i,/3i,a4) • E, 

B = -diag(/3i,/3i,-a4)"' -H. 

For an electromagnetic medium to be physically relevant, the Fresnel polynomial 
c^ijki should be hyperbolic polynomial [SWWIO]. This is a necessary con- 
dition for Maxwell's equations to form a predictive theory, that is, a necessary 
condition for Maxwell's equations to be solvable forward in time. The above 3D 
projections and the argument in [SWWIO] suggests that Metaclass II is hyperbolic 
for all /?! > while Metaclass IV is never hyperbolic for any Di g M. This is also 
supported by some numerical tests. 
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